Experimental analysis of decoherence in a CV bi— partite systems 
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Quantum properties are soon subject to decoherence once the quantum system interacts with the 
classical environment. In this paper we experimentally test how propagation losses, in a Gaussian 
channel, affect the bi-partite Gaussian entangled state generated by a sub-threshold type-II opti- 
cla parametric oscillator (OPO). Experimental results are discussed in terms of different quantum 
markers, as teleportation fidelity, quantum discord and mutual information, and continuous variable 
(CV) entanglement criteria. To analyse state properties we have retrieved the composite system 
covariance matrix by a single homodyne detector. We experimentally found that, even in presence 
of a strong decoherence, the generated state never disentangles and keeps breaking the quantum 
limit for the discord. This result proves that the class of CV entangled states discussed in this paper 
Oh would allow, in principle, to realize quantum teleportation over an infinitely long Gaussian channel. 
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"Quantumness" is very fragile and its surviving over long distances and times is subtly interconnected to the amount 
of noise coupling into the quantum system through its unavoidable interaction with the classical environment. Among 
the different quantum signatures, entanglement, a concept firstly introduced by Schrodinger [jj in response to the 
famous Einstein, Podolsky and Rosen paper in 1935 @, probably represents the most fascinating one. The conceptual 
consequences of entanglement are somehow subtle and deeply affect many aspects of modern physics as the interaction 
between reality and observer Q, the actual limit for defining metrology standards the limit to communication 
security @. 

The optical realization of such an intriguing quantum feature has gained, in the last decades, a central role for 
realising complex quantum communication protocols 6]. Among the different approaches the use of bright beams, 
OA ' i.e. optical modes carrying many photons, appears as the most robust carrier of entangled information [?J no matter 
if the latter will be coded into single-photon-like states (discrete quantum coding) or over the quadratures of the 
[**«- 1 carrier beams (continuous valued quantum coding) [8(. 

By interacting with the external world, a pure quantum state decoheres into a mixtur e [31 and its quantumness, 
if not completely death [l(|, need to be restored or re-enhanced by suitable protocols As a matter of fact, 

I [ decoherence limits both the attainable protocol fidelity and the type of accessible protocols. 
^ ■ The aim of this paper is to discuss and experimentally analyse the effects of the transmission over a lossy channel 
on the quantumness of bi-partite Gaussian continuous variable (CV) optical entangled states. While these effects on 
CV entangled states obtained by above threshold OPOs have been already investigated [l2[ , we focus our analysis on 
1 the states generated by a type-II sub-threshold OPO [13j. 

The birth of an entangled state happens whenever a genuine quantum correlation, i.e. without classical analogue, 
is set between two distinguishable quantum systems. This implies that the system wavefunction cannot factorize into 
the product of wavefunctions of the single sub-systems. While this definition is almost unique entanglement marks 
itself into different properties of the quantum state of the whole system and can be seen under different perspectives 
[Tij . For this reason, since the first experimental demonstration of entanglement [l5j several criteria and quantities 
have been introduced for analyzing entanglementness [lCil [l7| . Furthermore, the definition of an entangled state is 
conceptually connected to pure states ,18] . The properties of these states are described by wave functions and gedanken 
experiments can be though for testing different theoretical aspects but we all know that the states we can prepare 
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and manipulate in the laboratory are not pure but mixed. So that, experimental entanglement tests, especially in 
CV regime, are described in terms of density matrices rather than wavefunctions. Far from the ideal concept of pure 
states and systems that do not interact with the environment and/or observers, any quantum state undergoes to 
decoherence phenomena. In optics the most common process leading to decoherence is the phase insensitive loss of 
photons through diffusion and absorption mechanisms. This process is described by a Lindblad equation [19j for the 
evolution of the field operators that translates into a Master equation for the state density matrix. Eventually, if the 
state is Gaussian, i.e. its Wigner function is Gaussian, it admits a full description in term of its covariance matrix 

The investigated CV entangled state represents one of the possible quantum resources in CV teleportation protocols 
plj . The experimental data we present extend the analysis of Ref. [22| discussing the behaviour of the CV entangled 
system to the strong decoherence regime (up 99% of loss). Moreover, we discuss in details the relationship between 
the three different entanglement criteria used in the CV field linking them to the teleportation fidelity and quantum 
discord, two possib le q uantum signatures for evaluating the possibility of using this class of states in quantum commu- 
nication protocols |23l. I24I] . The discussion on the experimental results is preceded by a detailed theoretical foreword 
aimed at giving a physical insight the different criteria and quantities actually used for classifying CV entanglement. 
With this paper we prove that the state we are able to generate never disentangles. This rather counter-intuitive 
result has also been shown theoretically for the case of qubits states, showing a purity at the generation stage above 
a threshold value (25j . Moreover, the state we analyse can be used as the quantum resource for an effective quantum 
teleportation, of a coherent state, for any value of the transmission loss although the teleportation fidelity slowly 
reduce to its quantum limit. 

Optical entangled states can be obtained in non-linear processes such as parametric amplifiers that, dependin g on 
their operating regime, can prepare either single photon (2q | or continuous variable Gaussian entangled states (27l - l30j |. 
In the latter case the non-linear medium is allocated in a optical cavity and the OPOs, below the oscillation threshold 
and in a semiclassical approach, are described by bilinear Hamiltonians so realising the paradigm for Gaussian state 
generation 31] . In particular below threshold a single continuous wave OPO, generating squeezing in a fu lly dege nerate 
operation 32] , can give raise to a pair of bright CV entangled beams in the non-degenerate case [29l l30l . I331435} . Both 
the cases lead to states that represent robust resources for implementing different quantum communication tasks 36] . 

In this paper we experimental discuss the behaviour, under strong loss, of a bright bi-partite CV Gaussian entan- 
gled state outing a sub-threshold type-II OPO. In particular, we analyse, for this state, the behaviour of different 
quantumness and entanglement markers in order to discuss the limit at which this state can be transmitted before 
loosing its quantum ability of being employed in quantum communication protocols. More in details, we discuss: 
state purity, mutual information, quantum discord and three different entanglement criteria: Peres-Horodechi-Simon 
[37l HH, Duan [3^, and EPR-Reid In addition we relate these criteria to the teleportation fidelity F, i.e. the 
state ability of overcoming the quantum limit in a teleportation protocol for a coherent state. 

These states, already successfully generated in the laboratories [H, H3, HI] , are Gaussian, i.e. they are completely 
characterized by their covariance matrix, namely, by the first and second moments of their quadratures so that they 
can be easily characterized by a single homodyne detector [40] . 

The paper is organized as follows: first, in Sect. [IT] a general overview over the state generated by a sub-threshold 
OPO is given and the model for its evolution in a lossy Gaussian channel is discussed. In Sect. IIIII the different 
entanglement markers are introduced in terms of the state covariance matrix. The section includes a detailed discus- 
sion on the relationship between the different entanglement criteria and quantum markers for CV Gaussian states. 
Next, in Sect. IIVI the quantum markers' evolution in term of channel transmission is given. A brief description of 
the experimental apparatus, given in Sect. [Vj precedes the discussion, in Sect. IVI1 of our experimental findings. 
Eventually, in Sect. I VIII conclusions are drawn. 



II. GAUSSIAN CV BI-PARTITE SYSTEM 



A sub-threshold frequency degenerate type-II OPO generates a bi-partite state made of the two downconverted 
cross-polarized beams [41(. These states are Gaussian i.e. the Wigner function (W) of the state p, expressed in terms 
of field quadratures is Gaussian. In our case, we can assume that the average values of quadratures to be zero [l3| so 
that we can write W as 

exp [-iK T <T- 1 Kl 

where K = [Xt, Yi, X%, Y2) is the vector of the field quadratures (hereafter, whenever not necessary, we will omit 
the on operators so that: X = (a) + a ) /y/2 and Y = i ( a? — a) I v2) for mode 1 and 2 respectively. 
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A Gaussian state can be fully characterized by the matrix of the second order statistical moments of the field 
quadratures, i.e. the covariance matrix (CM) a, whose elements are defined by a^h = \ ({Kk, Kh}) ~ (Kk) (Kh) 
where {Kk, Kh] = K^K^ + K^K^. Moreover, for a bi-partite field state, a can be always given in the standard form 
by local symplectic operations 
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where a, (3 and 7 are 2x2 real matrices, representing respectively the self and mutual correlations matrices. The 
quantities, n, m, c\ and C2 are determined by the four local symplectic invariants 1\ = det(a) = n 2 , I2 = det(/3) = m , 
I3 = det(7) = C1C2, I4 = det(er) —Jnm — cf) (nm — c 2 ,). We note that the experimental matrices we have measured 
and that will be discussed in Sect. I VII are all in the above standard form. As a matter of fact, a sub-threshold type-II 
OPO, due to the symmetry of its Hamiltonian, can only produce states of this kind. So that, hereafter, whenever we 
refer to CM we will intend a covariance matrix in the standard form of Eq. ©. Moreover, at the time of their birth, 
states produced in an OPO show n = m and c\ = — C2, and the matrix is called symmetric and represents a bi-partite 
state where the energy is equally distributed between the two modes. 
<r is a bona fide CM i. e. it describes a physical state iff 

a + u> > , (3) 
where uj = f ^ J . The above condition can be written in terms of the four symplectic invariants 

h + I 2 + 2/3 < 4J 4 + i . (4) 

The above inequalities are equivalent to the Heisenberg uncertainty principle for a two-mode state and to ask covari- 
ance matrix positivity. 

The CM is also characterized by its symplectic eigenvalues 



d± 

Inequality (U) assumes a simple form in term of d- 



1 + I2 + 2/3 ± V Vi + h + 2/ 3 ) 2 - 4I 4 
\ 2 ' (5) 



We also note that a Gaussian pure state is a minimum uncertainty state and that the CM relative to a pure state 
necessary has det(er) = I4 = 1/16 (see Eq. (|11[) below) so that the case c\ — —c-i = c implies n = m, to ensure a bona 
fide CM, and for a pure fully symmetric state 



= V« 2 - 1/4 • (6) 



Moreover, for mixed fully symmetric states, c < \Jn 2 — 1/4 with the inequality saturated only by pure states. 
Hereafter we will indicate these fully symmetric states as diagonal. The meaning of this name will be clarified in Sect. 

inn 



A. Evolution in a lossy Gaussian channel 



Decoherence indicates the detrimental effect on a quantum system that stochastically interacts with the external 
world. The state of the system, initially pure, becomes mixed. In the density matrix language it translates into the fact 
that while a pure state is represented by an idempotent density matrix p (p 2 = p) this is not true for a decohered mixed 
state. The lossy transmission of an arbitrary optical quantum state between two sites is an irreversible decoherence 
process that can be described by using the open systems approach 19]. In this approach the environment is seen as 
a reservoir (thermal bath) made up of infinite modes at thermal equilibrium. The Kossakowski-Lindblad equation 
describes the time evolution in a noisy quantum channel of a multi-mode quantum state ps- This model is valid if 
the system S and the reservoir R satisfy the following general conditions: 
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• weak coupling (Born approximation) — the coupling between system and environment is so that the density 
matrix pp> of the environment is negligibly influenced by the interaction (the thermal bath state is stationary). 
This approximation allows to write the state psr (t) of the global system as psr (t) ~ ps (t) ® Pr; 

• Markovianity - there are not memory effects neither on the system and the reservoir. This approximation 
implies that the time scale r over which ps (t) changes appreciably under the influence of the bath is large 
compared to the time scale tr over which the bath forgets about its past, r 3> tr. 

• Secularity (rotating wave approximation) - the typical time scale ts of the intrinsic evolution of the system S is 
small compared to the relaxation time r. For an optical field it implies that the reservoir reacts to the average 
field and not to its instantaneous value. 

The evolution of an arbitrary two-mode Gaussian state in a noisy channel, at thermal equilibrium, can be translated 
into the formalism of the Wigner function obtaining the following Fokker-Planck equation for W (K) (see Eq.([T])) [H, 

d t W (K, t) = ~ (&.+^\ W (K) , (7) 

where T is the damping rate of the channel for both modes, while = (dx 1 ,dY 1 ,dx 27 dY 2 ) and V^- = d^dic = 

U X X + °Yi + U X 2 + U Y 2 ■ 

The evolution of Eq. ([7]) preserves the Gaussian character of the initial state and in terms of the CM a reads 

a(t)={l-e- rt )h + e- rt o-(0), (8) 

where <r(0) is the covariance matrix at t = and I is the 4x4 identity matrix (^1 is the vacuum state covariance 
matrix setting the standard quantum limit SQL). 

This form is in all equal to the effects of a fictitious beam-splitter (BS) that mimics the channel losses and couples 

into the system the vacuum quantum noise through its unfilled port. Being Uk (C) — e x P |C ( a l w fc ~ v k a kj | the 
SU(2) transformation induced by the BS on the fc-mode (fc = 1, 2, with Vk the modal operator for the vacuum) and 
T = e~ rt the power transmission of the beam splitter (tan£ = y/(l — T) /T), the above equation becomes: 

o T = (1 - T) h + Ta x ■ (9) 

In this form we can drop the temporal dependence and label the CM of the initial state as <jt=i = For states in 
the form of Eq. ([2|) the different elements evolve as 

m T = l+ T (m-~ 

Cl(2),T = Cx( 2 )T (10) 

For an infinite transmission channel T —> implies ot — > 5 1 i.e. the fully decohered CM represents a two-mode 
vacuum state. 



III. QUANTUM MARKERS AND ENTANGLEMENT 



The quantumness of a bi-partite CV state can be tested by two classes of markers. The first is intimately related to 
the state itself and includes inequalities that fix bounds for distinguishing among entangled and un-entangled states. 
The second has been translated into the quantum context from the classical information theory and is related to 
the amount of quantum information carried by the state. The latter includes quantitative measures such as mutual 
information, von Neumann entropy, quantum discord and teleportation fidelity. The former is made of entanglement 
criteria usually named by the authors that have theoretically found them. They are the PHS (Peres-Horodccki-Simon) 
[37], HH, Duan[39|. EPR-Reid [l6| criteria. Before going into details it is useful to start by introducing the state purity 
as an extra marker that tells how far the analysed state is from a pure one. 
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1. State purity 

Every quantum information protocol and every measurement scheme must take into account the fact that pure 
states are inevitably corrupted by the interaction with the environment and, therefore, in real experiments only 
mixed states are available. The purity p = Tr[p 2 ] (p < 1 for mixed state, = 1 for pure) may be expressed as a 
function of the CM © being 

(11) 



4 v /det [cr] 

Intimately related to the purity, it is possible to introduce the von Neumann entropy S (p). It is zero for pure states 
while strictly positive for mixed ones. In the case of the two-mode Gaussian state the entropy S (p) can be written 
in terms of its covariance matrix symplectic eigenvalues d± (see Eq. ([5])) as S (p) = f (d + ) + / [42| where 

/ (x) =(x + 1/2) log (x + 1/2) - (x - 1/2) log (x - 1/2) . (12) 

We note that both p and S (/?) are invariant under symplectic transformations. 



A. Entanglement criteria 

Although it has not been found a general solution to the problem of a quantitative measure of entanglement for 
mixed states, there exist necessary and sufficient conditions to asses whether a given state is entangled or not. These 
criteria provide a test for entanglement, and can be employed for studying the behaviour of entanglement in a quantum 
transmission channel. 



1. The Peres-Horodecki-Simon (PHS) criterion 

A bi-partite quantum state is separable iff its density operator can be written as a convex combination of the tensor 
product of density operators relative to the two different sub-systems [43[ 

p = ^2pjPji® Pj2, (13) 

3 

where y\ pj = 1 while pji i = 1, 2 are the density matrices of subsystems 1 and 2. By performing a partial transposition 
{i.e. transposition of the density matrix with respect to only one of the two Hilbert subspaces) p transforms into ppr- 
In the case p is a separable state, ppr should still represent a physical state so that ppr must be a non-negative 
density operator. On the contrary if ppr no more represents a physical state the system does not admit the form (|13|) . 
Thus, all separable states have a non-negative partially transposed density operator. From this consideration it is 
possible to deduce a necessary condition for separability or viceversa a sufficient condition for entanglement. In view 
of this the criterion is sometime referred to as the ppt criterion (positivity under partial transposition) [37ll38j. In the 
following we prefer to indicate it as the PHS criterion. It can be proven that it becomes a necessary and sufficient 
condition for Gaussian states [38| . From the CM point of view, partial transposition implies a sign flip for 1$ so that 
PHS criterion has a simple expression in terms of CM elements: a bi-partite Gaussian state is separable iff 

n 2 + m 2 + 2 \dc 2 \ - 4 (nm - c\) (nm - c 2 ) < - , (14) 

and it is entangled otherwise. We also note that the PHS criterion is invariant under symplectic transformations. 

PHS criterion relies on the possibility of describing independently the two subsystems. In phase space, transposition 
is defined as "time reversal" (or mirror reflection), which is given by a change of sign of out-of-phase quadrature 
(usually indicated as Y). Partial transposition is, therefore, a "local time reversal" which inverts the Y quadrature 
of only one subsystem. If any true quantum correlation is set between Yi and Y2 a sign flip on Y\ (or Y<i) will affect 
the sign of c\ (or C2) in Eq. ([2]) making ppt no more physical (see Eq. ((4])). 
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2. The Duan criterion 



For every bi-partitc state there exists a pair of EPR-like conjugate operators denned by 

u = \a\ x\ H — S?2 and v = \a\ p\ p 2 , (15) 

a a 

with a an arbitrary non-zero real number and = £8jj' (j,jf — 1, 2) and where subscript refers to the 1 (2) 

subsystem. Then, the birth of non-classical correlation between sub-systems 1 and 2 will lead to states for which 
the variance of EPR-like operators will reduce below the standard quantum limit (SQL). By calculating the total 
variance of such a pair of operators on p, a separable state of the form of Eq. (fT3|) . it can be proven [39| that 

<(A£) 2 ) p + ((A£) 2 ) p > a 2 + l, (16) 

setting a lower bound for separable states. Contrarily to the PHS criterion (see Sect. IIII A II) inequality (|16[) has 
been formulated only as a necessary condition for separability so that it is a sufficient condition for entanglement of 
a generic CV Gaussian state. 

As shown in [39| it becomes a necessary and sufficient condition for entangled CV Gaussian states. The sufficient 
and necessary condition can be expressed in terms of the covariance matrix elements iff the matrix itself is expressed 
in the form of Eq. (10) of Ref. [l| 

/ m a \ 

(17) 

with the matrix elements satisfying the constrains (11a) and (lib) of Ref. [39[ that, for the SQL equal to i, read 

r;i - 1/2 _ H2-1/2 
mi - 1/2 ~ m 2 - 1/2 ' 

|ci|-|c 2 | = VK - 1/2) (mi - 1/2) - y/fa - 1/2) (m 2 - 1/2) . (18) 
In this case the EPR operators pair of Eq. (|15l) are rewritten as 

u = a xi + sgn(ci) — x 2 and v = a pi - sgn(c 2 ) — p 2 , 
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where ao = y ^ = y ™ 2 2 _i/ 2 with the sufficient and necessary Duan criterion given by 

«o (ni + n 2 - 1) + — ^ 2 (Idl - |cb|) < . (19) 

a Q 

We note that, as proved in Ref. [39j. any Gaussian state can be transformed into the form (IT7|) by local linear unitary 
Bogoliubov operations, i.e. by acting independently on one or both the subsystems by applying local squeezing and/or 
rotations. 

In a more general fashion it is possible to write the sufficient but not necessary condition of Eq. for a CM in 
the standard form of Eq. ([2]) as 

(2n-l) a 2 + (2w 7 1) -2(c 1 -c 2 )<0, 



where a can be set to a 2 = y \^ to minimize the left hand side of the inequality: 



yj{2n - 1) (2m - 1) - (ci - c 2 ) < . (20) 



We note that, while for symmetric states (m = n) |a| = 1 and the EPR pair consists of two orthogonal field 
quadratures, this is not true, in general. 

The Duan criterion is strictly related to the Heisenberg principle for the single sub-system. If the state is separable 
the indeterminacy on a single operator is disjoint from the indeterminacy of the twin operator on the second sub- 
system; so that the total indeterminacy cannot violate the Heisenberg limit. This has nothing to do with conditional 
measurement and with the possibility of gaining information on one sub-system measuring the other. As we will see 
this approach leads to a stricter criterion: the so-called EPR " Reid" criterion. 
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3. The EPR "Reid" criterion 



A stronger bound can be found by analysing a bi-partite state under the shadow of conditional measurements. 
This concept descends directly from the original EPR gedanken experiment Q ■ For this reason it is usually indicated 
as the EPR criterion and was firstly introduced by Reid in 1989 [16j . in the very early days of quantum information. 
It describes the ability to deduce the expectation value of an observable on a sub-system by measuring the EPR 
companion observable on the second sub-system. The EPR "Reid" criterion sets, as the Duan one in its generic form 
([20| . only a sufficient condition for assessing entanglement. 

Mathematically this criterion can be deduced by calculating the conditional variance for an observable on sub- 
system 1 given the result of a measurement on sub-system 2 and comparing it with the standard quantum limit. 

For Gaussian states it can be written in terms of CM elements: 

\ nm ) \ nm J 4 

We note that, being based on conditional variances (and thus on conditional states) this last criterion is not symmetric 
under the exchange of the two sub-systems. This asymmetry reveals itself in a asymmetric sensitivity to loss pij . 
Then, the criterion itself can be recast if sub-system 1 is measured and the conditional variance on 2 is given 

\ nm J \ nm J 4 

The two definitions of the EPR criterion can make it ambiguous if one of the relations are not satisfied. This is not the 
case of balanced systems (m = n). Moreover, it can be proved that no pure state can asymmetrically violate the EPR 
criterion. At the same time this asymmetry can be useful for revealing eavesdropping in one-sided device-independent 
QKD schemes (45j . 

It is easy to see that the above two expressions for the EPR criterion are invariant for symplectic transformations 
like the PHS one (see Eq. ([T3|)). 

The EPR criterion is necessary and sufficient for assessing the "steering" form of nonlocality as recently discussed 
in Refs. 0II3]- This concept is very recently gaining much attention as a peculiar consequence of entaglement (46[. 
It is interesting to note that the concept of conditional variance, used for quantum steering, was used for assessing 
state preparation ability in quantum demolition experiments |47j. In that case steering happens between the probe 
and the signal systems. 



4- Witnesses 

All the above criteria (fT4")) , (|2"0|) . and (f2"T|) cannot be used other than as bounds. They are not suitable for measuring 
entanglement in a quantitative way. Quantifying entanglement for mixed states is a complicated issue still under 
discussion. There isn't, to date, a single, universal measure that quantifies the entanglement for a mixed state. 

In discussing the experimental results in Sect. I VII we will make use of three different witnesses, one each for the 
above criteria, in order to evaluate their behaviours under lossy transmission: 

wphs = 4 (nm — cf) (nm — c 2 .) + - — (n 2 + m 2 ) — 2 \c\Cz\ , 



wduan = m \n - - m - - - (ci - c 2 




WEPR. = ^1-^1-^ --. (22) 



The three w's don't satisfy the requirements for being a measure of entanglement. For example they don't verify 
the basic axiom stating that a good measure should be equal to for any separable state (48|. 

In summary, for a Gaussian bi-partite state the three criteria (see Eqs. ([T4|) . (|20]l. and ([2T]) ) reduce to 

wphs < 

p is entangled ^ J wduan < . (23) 

WEPR < 
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However, once the state p is entangled wpns, wduan and wepr provide suitable markers for evaluating how far 
the state is from being separable, somehow measuring the robustness of the entanglement. 

Eventually we note that for diagonal fully symmetric states (n = m and c\ = — C2 = c in Eq. ([2])) wphs, wduan 
and wepr read 



wphs = 
wduan = 

wepr = 




and the two bounds (c > n — 1/2) for wpns and wduan coincide while the bound for wepr. is c > y n { n ~~2~) so 
that the EPR criterion is stricter than the PHS and Duan ones for any allowed value of n. 

B. Quantifying entanglement 

The quality of an entangled state can be, indeed, evaluated, operatively, by looking at the effectiveness such a state 
could, in principle, give in quantum communication protocols. To this end we use J 7 , the fidelity of CV teleportation 
of a coherent state, as a benchmark. In fact, J 7 for a Gaussian resource depends only on the entanglement qualit y of 
the resource itself. Using as resource a Gaussian two-mode state described by the CM in Eq. (J2|), T is given by [49| 

T = (1 + m + n - 2ci)~2 (i + m + n + 2c 2 )~^ . (24) 

It is possible to show that a teleportation protocol of a coherent state, fully based on classical strategies, provides 
T <l/2. Then, T >l/2 implies that the resource is entangled. So that, the value of F becomes an indicator of the 
quality of the entanglement. We note that for diagonal states (n = m and c\ = —C2 = c) J- = (1 + 2n — 2c) _1 = 
(2 + wduan) ■ In this particular case, perfect fidelity would be obtained for wduan = —2 that would imply an 
unphysical CM so that Gaussian resources, as the state produced by OPOs, cannot guarantee perfect teleportation 
[5(J |. Moreover, T > 1/2 gives c > n — 1/2 and coincides with both the Duan and PHS bounds (see the end of the 
previous subsection) for such diagonal states. We note that, similarly to the Duan (Eq. PU])') but contrarily to the 
PHS and EPR criteria (Eqs. (Q~4|) and ([21]) ) T is not invariant under symplectic transformations. 

As already mentioned, only the PHS criterion can be written for a generic CM as a sufficient and necessary 
condition while the Duan and the EPR ones set only sufficient bounds. In order to discuss the relations between 
the different bounds and the teleportation fidelity, as an entanglement marker, it is possible to draw a region plot, 
similar to the one drawn in Ref. [5l[ for a different case. In Fig. [TJ we have visualized the different bounds set by 
the three entanglement witnesses of Eq. (|2"2"]0 and the region for which T > 1/2 (see Eq. (|2"4")l ). The plot has been 
computed considering a CM in the form of Eq. ([2j with m — n (balanced system). The axes report the value of 
Si (below) and c 2 (left), the two correlation terms of the covariance matrix are normalized to cm ax = \/n 2 — 1/4 
so that ci = — C2 = 1 represents a pure maximally entangled state (i.e. the state showing the maximum quantum 
correlation for a given total energy of the system). We note that fully symmetric states, that we indicated as diagonal 
states, lay on the plot diagonal (top-left to bottom-right). These states, besides their symmetry, can be obtained at 
the end of a lossy propagation, described by Eq. ©, of an initially pure state. On the contrary, states outside the 
diagonal, having c\ ^ — £2, cannot be obtained by propagating pure states. As a matter of fact, CV entangled states 
produced by type-II OPO, show, in view of the symmetry in the interaction Hamiltonian, ci = —£2. This is not true 
for CV entangled states obtained by mixing the outputs of two independent type-I OPOs. In such a case the two 
fields have disjoint Hamiltonians and the symmetry is broken. 

The light gray (labelled with (VI)) areas indicate un-physical states i.e. CMs violating inequality (|4]). 

Diagonal states satisfy the conditions (see Eq. (|18p) for which the Duan criterion becomes also necessary so that the 
coincidence between the Duan and the PHS bounds, along the diagonal, is not a surprise. Being both necessary and 
sufficient they coincide. For these diagonal states entanglement (seen as non-separability property) implies J- > 1/2. 
So that entanglement is a pre-requisite for using the state as a resource for the teleportation of a coherent state. 

There are two interesting regions in the plot that deserve some comments. The first one, encompassing areas labelled 
as (777) and (IV) in the plot (lights green and blue), represents states that violate the PHS bound (wpns < 0) while 
they do not the Duan one (wduan > 0). This apparent ambiguity can be solved noting that the CMs relative to 
these regions do not respect the condition (|18j) so that the non-negativity of wduan does not imply a separability 
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Cl 

FIG. 1: Region plot [color online) of the different entaglement witnesses of Eqs. \22\ and teleportation fidelity (Eq. (|24l) ~) as 
an entanglement marker. The light gray (labelled with (VI)) areas indicate un-physical CMs (i.e. violating inequality ©). 
The different criteria show different regions of entanglement (see text for details). 

of the state. On the other hand for such states wpns < implies that they are effectively entangled and that their 
density matrix p cannot factorize into a convex combination of the tensor product of density operators relative to 
the two different sub-systems. It is possible to see that if these CMs are transformed by local squeezing operations, 
as outlined in Ref. [39], into a form that fulfills conditions (Tl8l) . the transformed states show wduan < 0. We have 
numerically done a few tests on such odd matrices, verifying that once taken into that form, the states violate the 
Duan bound (fT9")) . We note that, being the latter written in a more general form, it is more useful from the practical 
point of view. 

The second interesting region, labelled with (IV) (light green) in Fig. [TJ represents states that, although entangled, 
cannot be used for teleporting coherent states. CMs lying inside this area will not give T > 1/2. It is interesting to 
note that such states fall also inside the region for which the Duan criterion (I2U1) is not fulfilled. As above mentioned, 
once the relative CMs are transformed into the form (TP7|) °v local squeezing the transformed state will fulfill the 
Duan criterion in the form (|19p so that, in this new scenario, the system will be entangled. At the same time, if this 
novel state is used as a quantum resource for teleportation of a coherent state it will give J- > 1/2 [4!| so that local 
squeezings unveil entanglement. The initial state lying in this area is entangled for PHS, being wphs < 0, but, from 
the point of view of teleportation, entanglement manifests itself in an useless way. This entanglement can be made 
useful by locally transforming the two subsystems. 

We can see that the EPR criterion (region (/), light brown) offers a more restrictive condition with respect to the 
other two criteria even for diagonal states. 

Region (II) (salmon) represents the bound fixed by the Duan criterion as a sufficient but not necessary condition. 
While region (V) (white) represent separable states. 

1. Mutual information and quantum discord 

Recent studies have shown that some particular classes of separable correlated states, traditionally considered 
classical, show quantum features useful for application in quantum technology p52l |. 

In statistical theory any correlation between two random variables A and B can be measured by their mutual 
information defined by two equivalent expressions: 



I(A;B) = H(A)+H(B)-H(A 7 B) 

= H (A) — H (A\B) = H (B) — H (B\A) 



(25) 
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where H (X) is the Shannon entropy (the classical analogue of the S (p), von Neumann entropy defined in Sect. IIII 1[) 
associated to the statistical distribution of the variable X ((X, Y) indicating the joint probability distribution and 
(A|Y") the conditional distribution of Y given the value of X). These two definitions, translated into the quantum 
language by substituting H with S (p), the von Neumann entropy, do not coincide anymore [53| . So, while the 
translation of the first one into the quantum language is straightforward and univocal, this is not true for the second 
one. 

The first of the two definitions (|25[) . is unambiguously referred to as the quantum mutual information [54| between 
state 1 and 2 (p representing the state of the bi-partite system as a whole) 

l(p) = S(pi) + S( P2 )-S(p) . (26) 

where pi(2) = Tr 2 (\) [p] are the partial traces. It relies only on the quantum state of the whole system compared 
with the single sub-systems' states. In a bi-partite system, described by a density matrix p, T(p) quantifies the total 
correlation between the subsystems pi and p^. It is > for entangled states, while it is strictly = for separable 
systems. It can be written in terms of the CM invariants (see Sect. [TT]) and symplectic eigenvalues (see Eq. ([5])) 

X (a) = f (n) + f (to) - / (d+) - / (d_) , (27) 

with / (x) given in Eq. (Ti"2l . 

The second definition of Eq. (|25p. translated into the quantum world, necessarily involves the conditional state of 
a subsystem after a measurement performed on the other one. So that, the symmetry between the two subsystems is 
broken. Since the conditional entropy H {A\B) requires to specify the state of B given the state of A, its definition, 
in quantum theory is ambiguous until the to-be-measured observables on A are selected so that the conditional state 
of B can be defined. 

This discrepancy has lead to the concept of quantum discord T> (p) [55j . A non zero V (p) signals the presence 
of quantum features in the correlation between the two sub-systems notwithstanding their separability or entangled 
nature. So that T> (p) is a measure of genuine quantum correlation 

V(p)=l(p)-C(p), 

where C (p) is the amount of genuinely classical correlation. 



For a Gaussian state described by the CM ([2]), T> becomes 53 1 



V (a) =/(to)-/ ( d+ ) - f (d_) + / ( n + 2 1 7 2 + m 2Cl ° 2 ) • (28) 

We note that, as for the EPR criterion, in quantum discord there is an asymmetry in the exchange of the two 
sub-systems. Again this is due to the use of the concept of conditional states. 

IV. QUANTUM MARKERS EVOLUTION 

As shown the CM of a bi-partite state undergoing to a lossy transmission evolves as Eq. Consequently the 
different quantum markers evolve as: 



Pt 



4-y/det [cr T ] 



1 + (tot + nr) 2 + 2 (c2,t — ci,r) (1 + «t + nr) + 2 (tot + nr — 2ci,tC2,t) 



w phs,t = u t + m T + 2 |ci,tC2,t I — 4 [n^vrij- — c( T ) (jiTrriT — c^t) 
wduan.t = \J (2n T - 1) (2m T - 1) - (ci t - c 2 ,t) 



wepr,t 




c. 



2 

2,'i' 



\ 1 



riTTriT I \ riTrriT 



4 



1 T = I (n T ) + f (to t ) - / (d+, T ) - / (d-, T ) 

v T = /(tot)-/ ( d+ , T) - f (d_, T) + / ( " T+2 ;nC 2ciC2 ) (29) 



11 



where the subscript t indicates the quantity undergone to a lossy transmission described by Eq.©. 

We note that the vacuum state obtained for T = is a pure one i.e. [1q = 1. Moreover, in the ideal case (no loss), 
the OPO would generate a pure state as well. Being /Mr < A*o,i) for T ^ 0,1, the evolution of /it is not monotonic. 
The purity of the composite system cannot be considered a general entanglement marker (56[. As a matter of fact, 
any pair of physical systems in a pure state have \i — 1 even if they are disentangled. On the other hand, in our 
specific case, having a precise hypothesis on the ideal state (a pure twin-beam diagonal state) outing the OPO crystal 
allows us to consider [i as a measure of the total decoherence that has affected the state. 

It is easy to see that Tt, wphs,t, and u>duan,t describe properties very robust under decoherence. Once T > 1/2, 
wphs < an d wduan < for T = 1 they will keep breaking the respective bounds for every level of loss. Both 
mutual information and quantum discord show, with respect to loss, the same feature even if decoherence affects their 
amount. On the contrary a state that show EPR-like correlation (wepr < 0) for T = 1 will not keep this property 
along the propagation so that some particular protocol based on this property cannot be reproduced. Under a total 
loss greater than 50% any state looses this correlation property. 



The transmission over a Gaussian channel is described by Eq. ((TJ- As already mentioned this evolution is in all 
equivalent to a fixed amount of loss introduced by a fictitious beam-splitter. The actual experimental apparatus is 
made of the CV entangled state source, a variable attenuator (mimicking the BS), and a state characterization stage. 
A block diagram of the experimental setup is presented in Fig. [2] 



FIG. 2: (color online) Block diagram of the experimental setup. The details on the OPO source are given in Ref. 57], while 
the characterization stage, based on a single homodyne detector, is fully described in Ref. [40| ]. 



The state source relies on a CW internally frequency doubled Nd:YAG laser (Innolight Diabolo) whose outputs 
@532nm and @1064nm are respectively used as the pump for a non degenerate OPO and the local oscillator (LO) 
for the homodyne detector. The OPO is set to work below the oscillation threshold and it provides at its output two 
entangled thermal states (the signal a and the idler b). 

The OPO is based on an a-cut periodically poled KTP non linear crystal (PPKTP, Raicol Crystals Ltd. on custom 
design) which allows implementing a type II phase matching with frequency degenerate and cross polarized signal and 
idler beams, for a crystal temperature of ~ 53°C. The transmissivity of the cavity output mirror, T outl is chosen in 
order to guarantee, together with crystal losses (re) and other losses mechanisms (T;„), an output coupling parameter 
Vout — Tout/(Tin + re) @1064 nm of ~ 0.73, corresponding to an experimental linewidth of 15 MHz @1064 nm. In 
order to obtain a low oscillation threshold, OPO cavity geometry is set to warrant simultaneous resonance on the 
pump, the signal and the idler {57]. Measured oscillation threshold is around 50 mW; during measurement runs the 
system has been operated at «60% of the threshold power to avoid unwanted non-Gaussian effects (5|| . 

The two beams outing the OPO are transmitted through a filter of variable optical density that mimics the BS. 
The loss level introduced by the filter is polarization independent and can be tuned from a few percent up to more 
than 99%. 



V. THE EXPERIMENT 




A. CV entangled state source 
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B. Characterization stage 



The signal and idler modes are then sent to the covariance matrix measurement set-up: this consists in a preliminary 
polarization system, that allows choosing the beam to be detected and a standard homodyne detector. The polarization 
system is made of an half-wave plate (A/2) followed by a polarizing beam splitter (PBS); the different wave-plate 
orientations allow choosing the beam to be transmitted by the PBS: the signal (a), the idler (b) or their combinations 
c = -±=(a + b) or d = -j^{a — b). Two other combinations, e = -^(ia + b) and / = -^=(ia — b), can be obtained by 

inserting before the PBS an additional quarter wave plate (A/4) [40|. Acquisition times are considerably short thank 
to pc-driven mechanical actuators that allow setting the A/2 and A/4 positions in a fast and well calibrated manner. 

Once a beam is selected, it goes to an homodyne detector put downstream the PBS. This exploits, as local oscillator, 
the laser output @1064 nm, previously filtered and adjusted to match the geometrical properties of the OPO output: 
a typical interferometer visibility is 0.98. The LO phase 9 is spanned to obtain a 2ir variation in 200 ms. In order to 
avoid low frequency noise the homodyne current is demodulated at £1=3 MHz and low-pass filtered (i?=300 KHz). 
Then it is sampled by a PCI acquisition board obtaining 10 6 pts/run, with 14-bit resolution. The electronic noise floor 
is 16 dBm below the shot noise level, corresponding to SNR w40. Data are analysed by a ©Mathematica routine 
that extract from the six homodyne traces the relevant quadrature variances necessary for reconstructing the whole 
covariance matrices [40l |. 



VI. EXPERIMENTAL RESULTS 



We have performed different sets of measurement under lossy transmission in order to investigate different loss 
regimes. Each experimental CM comes from seven homodyne traces: the shot noise calibration trace, obtained by 
obscuring the OPO output, six traces each for one of the six modes required for a full state characterization. Contrarily 
to other previous experiments, where quantum tomographic routine were used in order to retrieve experimental CMs 
[l3[, we have evaluated CMs by a simpler ©Mathematica routine that calculates relevant second order moments of 
homodyne distributions in a faster way without enhancing the experimental indeterminacy on the CM elements. We 
have tested on a few CMs that this procedure gives results in all compatible with quantum tomography. We have 
also checked, with the standard procedure outlined in 59], that the states under scrutiny were effectively Gaussian. 

Once a CM is obtained the different entanglement witnesses (ifpirs, wduan, and wepr Eqs. (|22p ). state purity 
(/i Eq. (JTTJ)), teleportation fidelity (F Eq. ([M])). quantum discord (V Eq. ([2"5])). and mutual information (1 Eq. (|2"o]l ) 
are calculated. Then, the overall decoherence, i.e. the total level of loss that includes OPO cavity escape efficiency, 



propagation loss, filter absorption, homodyne efficiency, is assigned as a label to the measurement [22j. This expected 
level of decoherence is then compared to the theoretical one obtained by inverting Eq. ^ and solving for T under 
the condition det(eri) = 1/16; thus requiring that o\ represents a pure state (see the discussion at the end of Sect. 
HI]l . We have verified that, even if experimental CMs do not reproduce exactly diagonal states, all of them respect 
the Duan conditions (|18[) within experimental indeterminacies. So that for the analysed matrices the Duan witness 
wduan represents a sufficient and necessary condition for entanglement. 

In all the reported plots we have considered the less decohered datum (obtained for T = 0.63) as a reference so 
that all the reported theoretical curves are obtained imposing that Eqs. (fTUl) and (|2"9"|) evaluted for T — 0.63 give the 
measured values. 

The total losses we have measured span the interval 37 — 99% (0.01 < T < 0.63). We note that T — 0.63, in absence 
of extra loss and having a cavity escape efficiency of « 0.73, implies an overall state detection efficiency of ~ 0.86 in 
agreement with an homodyne visibility of 0.98 ± 0.01, a photodiode (nominal) efficiency of 0.90 ± 0.01 and residual 
transmission loss between the OPO output mirror and the detector surface of 0.01 ± 0.01. 

In Fig. [3] we report the behaviour vs. T of the averaged correlation term ((|ci,t| + |c2,t|) /2 see Eq. (flU)) ). As 
expected the correlation between the two sub-systems degrades linearly with the total loss (T — > 0). The expected 
behaviour (full dark orange line), obtained by considering the less absorbed CM (T = 0.63) as a reference, follows 
quite well the reported data. Actually, data refer to acquisition taken on different days so that, the scattering of the 
point around that line is more due to source long-term dynamics then to actual deviation from the Lindblad model. 
At the same time the fact that the points are reasonably close to that line proves that the long term stability of the 
source can be considered quite good. 

As already mentioned, wpns an d wduan describe a physical property of the state that is strong under decoherence 
as proved for lower loss (below 90%) in Ref. [22| ■ They are symptoms of un-separability, in the sense that the system 
state cannot be described by a density matrix in the form of Eq. (fT"3|) [43]. This can be seen in Figs. [4] and [5] where 
wphs and wduan are plotted vs. T. We have also enlarged, in the insets the region for strong loss (T < 0.15) 
to prove that, even if the analysed state is very close to a two mode vacuum (the total average number of photon 
((n + m — 1) /2) reduces to 0.02 ± 0.01) it is still experimentally possible to prove that the state is non-separable. It 
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FIG. 3: (color online) behaviour of the averaged correlation term |ci,t| + |c2,t| /2 in Eq. ([9]). As expected the correlation reduces 
linearly with T. The full [dark orange) line represents the linear behaviour calculated starting from the first experimental point 
we have measured (T = 0.63). Error bars are smaller than data points and amount to ±0.01. 
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FIG. 4: (color online) wphs vs. T. The full (dark orange) line represents the expected behaviour calculated by the third 
of Eqs. Q29p setting the first experimental point at T = 0.63 as the intial datum. Error bars, obtained by propagating the 
experimental indeterminacies in Eq. (|22a l. range between 10 -4 and 0.1. In the inset we report the high loss regime (T < 0.15) 
for better enlightning the un-separability, as witnessed by wphs, even in presence of strong decoherence. 
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FIG. 5: (color online) Wduan vs. T. The full (dark orange) line represents the expected behaviour calculated by the fourth 
of Eqs. (|29(l setting the first experimental point at T = 0.63 as the intial datum. Error bars, obtained by propagating the 
experimental indeterminacies in Eq. (|22b ). range between 0.01 and 0.06. In the inset we report the plot for T < 0.15 in order 
to better visualize the persistence of entanglement, as witnessed by wduan, even in presence of strong decoherence. 
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has to be noted that, while, for T — > 0, wpns approaches its classical limit non-linearly (see the third of Eqs. (|29]l ). 
wduan (see the fourth of Eqs. (|29p ) is linear. Thus, in the very high loss regime it becomes more reliable to assess 
entanglement using the latter than the former. 



0.1 
0.0 

wepr 

-0.1 
-0.2 



0.0 0.2 0.4 0.6 0.8 1.0 

T 

FIG. 6: (color online) wepr vs. T. The full (dark orange) line represents the expected behaviour calculated by the fifth 
of Eqs. (|29[) setting the first experimental point at T = 0.63 as the intial datum. Error bars, obtained by propagating the 
experimental indeterminacies in Eq. (|22b ). range between 2 x 10~ 4 and 0.02. They are considerably larger for point at low 
losses. As expected for total losses larger than 0.5 wepr > and the state does not show EPR correlation. 

Wepr < indicates that the state exhibits EPR-like correlation so that it is possible to gain information on the 
expectation value of one observable on one sub-system with a precision better than the standard quantum limit once 
the EPR companion is measured on the other sub-system. This feature is by far the most fragile under decoherence: 
for T < 0.5 no state can keep this quantum feature. This can be understood from the fact that loss, a stochastic 
process, affects directly the degree of correlation between the two sub-systems while the system representation (i.e. its 
un-separability) is only smoothed by this process. It is relevant to note that T = 0.5 also corresponds to the minimum 
state purity \i. So that, loosing the EPR character coincides with the maximum mixedness for the state during its 
propagation. In Fig. [5] we report the experimental behaviour found for wepr for our system. Measured CMs for 
T < 0.5 all show wepr > 0. A positive wepr. indicates that, for these states, any attempt to gain information on 
one sub-system by measuring the other would result less precise than the standard quantum limit. 
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FIG. 7: (color online) T vs. T. The full (dark orange) line represents the expected behaviour calculated by the second of Eqs. 
(|29|l setting the first experimental point at T — 0.63 as the intial datum. Error bars, obtained by propagating the experimental 
indeterminacies in Eq. (|24|) . range between 10 -4 and 0.01. In the inset we report the plot for T < 0.15 in order to underline 
the persistence of a quantum teleportation regime even in presence of strong decoherence (high loss) thus proving that even 
for, in principle, infinite distance this class of states would allow to perform the teleportation of a coherent state with a fidelity 
above 1/2. 

An important signature for an entangled CV state is its ability of acting as a quantum resource in the CV telepor- 
tation protocol for coherent state. In Eq. (|24l) we have expressed the fidelity J 7 as a function of the CM elements. 
J 7 , as wpns a n d wduan, represents a robust signature of quantum properties for the state undergoing to a lossy 
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transmission. In particular, in Fig. [7j we see that even in the high loss regime, J- remains above the classical limit 
of 0.5 (see the inset for greater details). Thus proving that CV entangled state, as the one produced by our source, 
could be used as resource for realising teleportation protocol of coherent state, in principle, at an infinite distance. 

Eventually we have retrieved, from our CMs, the value for the quantum mutual information 1(a) (Eq. (|27p) and 
quantum discord V (a) (Eq. p5|)). 
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FIG. 8: (color online) T and T> vs.T. The full (dark orange) and dashed (blue) lines represent the expected behaviours 
calculated by the sixth and seventh of Eqs. (|29p setting the first experimental point at T = 0.63 as the intial datum. Error 
bars, obtained by propagating the experimental indeterminacies in Eqs. (|27[) and (|28[) , respectively, range between 3 x 1CP 3 
and 0.02 for 1 and 10" 4 and 0.03 for V. In the inset we report the V data for T < 0.15 in order to underline the persistence 
of true quantum correlation even in presence of strong decoherence (high loss) . Note that the data for I scatter more from the 
expected behaviour may be signalling extra classical correlations. 

In Fig. [8] we report the experimental data together with the expected behaviours, as usually calculated considering 
the less decohered matrix as a reference, for X and T> vs. T. As it can be seen the quantum discord follows very 
well its theoretical line while quantum mutual information is a little more scattered around it. Moreover, our data 
prove that even in presence of strong decoherence, it is possible to evaluate that T> keeps > 0, within the experimental 
indeterminacies, all the way down to an highly absorbed state. We note that Gaussian quantum discord is attracting, 
very recently, a lot of experimental interest O HTj, [6l[. In particular, in Ref. [24j . the authors give an operational 
significance to quantum discord as the possibility of encoding quantum information in separable states. In Ref. (60} 
the optimal strategy for evaluating T> (a) in homodyne measurement is presented. It is interesting to compare our 
experimental plot with the one reported in Ref. (6i| where the authors analyse the quantum discord under the lossy 
transmission of one of the two sub-systems. We note that in their case the scattering of the experimental points 
around the theoretical curve is almost equivalent for I and T> while in our case there is a clear difference. 

VII. CONCLUSIONS 

Gaussian bi-partite states are one of the most renown resources for implementing CV quantum communication 
protocols such as CV teleportation of coherent states. In this paper we have experimentally analysed how decoherence 
affects different entanglement criteria and quantum markers for a CV bi-partite state outing a sub-threshold type-II 
OPO. The decoherence is experimentally introduced by transmistting the quantum state through a variable attenuator. 
Before illustrating our experimental results we have discussed in details the relationship between the three different 
entanglement criteria used in the CV framework and linked them to the teleportation fidelity and quantum discord. 
The latters represent two possible quantum signatures for evaluating the ability of this class of states in quantum 
communication protocols. 

On one hand, our findings prove that the Lindblad approach for describing lossy transmission is valid all the way 
down to strongly decohered states. On the other hand, with this paper, we prove that the particular class of states 
we have analysed keeps, within the experimental indeterminacies, its main quantum signatures, i.e. the possibility of 
realizing quantum teleportation of coherent states with a fidelity above 0.5 and a quantum discord above for a total 
loss of ~ 99%. This proves that the class of CV entangled states, we analysed, would allow, in principle, to realize 
quantum teleportation over an infinitely long Gaussian channel. 
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In analysing how quantum discord (see Fig. [8]) and quantum mutual information behave under decoherence we 
interestingly found that the scattering of the points around the theoretical curve is significantly more evident for 
the quantum mutual information may be signalling that a key role, in our case, is played by unexpected classical 
correlations. This point will be subject of further theoretical and experimental investigation. 

The authors thank S. De Siena and F. Illuminati for useful suggestions and discussions on the theoretical aspects. 
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